Olj 


UDC-n-78>U3  hM  bMB  rwrlMMd  aad  U •ptmrrmd  for  poblicatlon. 


'ikast  SCHEBBI 
Projset  Wngln— r 


APPIOVBDt 


dS 


JOSEPH  J.  MABBSnr 

Chl«f,  lAllablllty  and  Co^Mtibllltj  Dlvlsi^ 


POE  THE  OdtIARDER: 


JOHN  P.  HDSsM'id'' 

Acting  Chief,  Plans  Office 


TO  rmpen  bm  »aaa  rmrUmmJ  by  tha  BADC  Infozmtlon  Office  (01)  and  is 
Mlaaa^la  to  latic^  Tarliical  Znfonation  Sarrica  (NTIS).  At  NTIS  it 
•ill  ba  ralaaaabla  to  E— ral  poblic,  including  foraigB  nations* 


t*  b.  rwv.d  fro.  th.  MDC 

aaintaining  a currant  MUing  list.  •■■lat  ua  in 


Bo  not  ratum  this  copy,  latain  or  dastroy. 


UNCLASSIFIED 


iSSiriCATION  or  ThK  PAOC  (Whu  Dl»  Knltn4) 


y REPORT  DOCUMENTATION  PAGE 


BEFORE  COMPLETTOG  FORII 
RKClPICNT^t  CATAtOO  HUMFCR 


00 VT  ACCKUION  NO. 


lY  SJRUCTUSE  OP  DUCTED  MULTIPATH  FOR  AN 
IS  LINK- 


1.  AUTHOIV*; 


« rcnromtiNO  onoanization  name  ano  adoness 
George  Ueshington  University 

Dept  of  Electricel  Engineering  & Computer  Science 


Wash  DC  20052 


II.  contnollino  orricE  name  ano  aodneei 


Rone  Air  Development  Center  (RBC) 

Grifflss  APB  HY  13AA1 ! 

tnSNTYoRmSTTSISFvNAMrSAOONEiiTlTdlfJIflnlfroiiiCoJilrolTftiioWeiJ 


IS.  SECUNITV  CL 


Rome  Air  Development  Center  (RBC) 
Grifflss  APB  NY  13441 


IS*.  OECLASSIFICATION/DOWNGRAOING 
SCHEDULE 


Approved  for  public  release;  distribution 


IS.  sufflementany  notes 


RADC  Project  Engineer:  Jacob  (NMI)  Scherer 


IS.  KEY  WONOS  fConlinu*  on  roYoro*  oltfo  If  noeooaarr  ontf  lOonlIfy  ky  btocA  nuaibor; 


Propagation 

Llne-of-slght 

ConBBunlcatlons 


20.  ASSTMACT  (Conttnu*  on  roooeoo  tidm  U noeoooory  md  IdontHy  by  block  numbof) 


This  report  describes  the  geometric  optic  analysis  of  a duct  model  for  llne- 
of-slght  (LOS)  fading.  The  atmospheric  refractive  Index  Is  modeled  by 
three  linear  segments  so  as  to  form  an  elevated  duct.  The  earth's  curvature 
is  taken  into  account  by  using  the  modified  refractive  Index.  Simple 
approximate  formula  are  developed  for  the  number  of  geometrlc-optlc  rays 
received  for  a given  link  length.  In  addition,  the  angle  of  departure  from 
the  transmitter  of  these  rays  Is  calculated.  The  angles  enable  one  to  . 


aaCuMTV  CLASSIFICATION  OF  THIS  FAOE  (Vliwi  OMa  SntarsO 


UNCLASSIFIED 


rMdily  calculata  th«  uplltuda  and  phaaa  of  contributing  raya.  To  gain 
Inalght  Into  tha  atructura  of  tha  ray  faalllaa,  tha  trajactorlaa  of  trappad 
raya  ara  plottad.  Wo  hava  found  that  tuo  cuapa  axlat,  aach  cuap  having  tuo 
cauatlca  aaaoclatod  with  It. ^ 


(tCuniTV  CkAUiriCATION  OA  TMIt  PAOCfVkM  0«A 


TABLE  OF  CONTENTS 


PAGE 


I  IntrcxJuction  1 


II  Channel  Modeling  3 


III  Multipath  Ray  Structure  for  a LOS  Link  in  the 


Presence  of  a Duct.  8 

A.  Problem  Formulation-Modified  8 

Refractive  Index. 

B.  Geometric-Optics  11 

C.  Direct  Ray  Calculations  15 

D.  Multipath  Rays  18 

E.  Ray  Families  24 

IV  Conclusions  and  Recommendations  32 

References  34 


1 


PREFACE 

nils  effort  was  conducted  by  George  Washington  University  under  the 
sponsorship  of  the  Rome  Air  Development  Center  Post-Doctoral  Program  for 
Rome  All  Development  Center.  Mr.  Joseph  Mensch'of  DCA  was  the  task  project 
engineer  and  provided  overall  technical  direction  and  guidance. 

The  RADC  Post-Doctoral  Program  Is  a cooperative  venture  between  RADC 
and  some  slxty-flve  universities  eligible  to  participate  in  the  program. 
Syracuse  University  (Department  of  Electrical  Engineering),  Purdue  Univer- 
sity (School  of  Electrical  Engineering) , Georgia  Institute  of  Technology 
(School  of  Electrical  Engineering),  and  State  University  of  New  York  at 
Buffalo  (Department  of  Electrical  Engineering)  act  as  prime  contractor 
schools  with'  other  schools  participating  via  sub-contracts  with  prime 
schools.  The  U.S.  Air  Force  Academy  (Department  of  Electrical  Engineering), 
Air  Force  Institute  of  Technology  (Department  of  Electrical  Engineering) , 


and  the  Naval  Post  Graduate  School  (Department  of  Electrical  Engineering) 
also  participate  in  the  program. 

The  Post-Doctoral  Program  provides  an  opportunity  for  faculty  at 


participating  universities  to  spend  up  to  one  year  full  time  on  explora- 
tory development  and  problem-solving  efforts  with  the  post-doctorals 
splitting  their  time  between  the  customer  location  and  their  educational 
Institutions.  The  program  Is  totally  customer-funded  with  current  pro- 
jects being  undertaken  for  Rome  Air  Development  Center  (RADC),  Space  and 


Missile  Systems  Organization  (SAMSO),  Aeronautical  System  Division  (ASD) , 
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ill 


Electronics  Systems  Division  (ESD),  Air  Force  Avionics  Laboratory  (AFAL) , 
Foreign  Technology  Division  (FTD),  Air  Force  Weapons  Laboratory  (AFWL), 
Armament  Development  and  Test  Center  (ADTC) , Air  Force  Conmunlcations 
Service  (AFCS),  Aerospace  Defense  Command  (ADC),  HQ  USAF,  Defense  Com- 
munications Agency  (DCA),  Navy,  Army,  Aerospace  Medical  Division  (AMD), 
and  Federal  Aviation  Administration  (FAA) . 

Further  information  about  the  RAOC  Post-Doctoral  Program  can  be 
obtained  from  Mr.  Jacob  Scherer,  RADC/RBC,  Grlfflos  AFB,  NY  13441, 
telephone  Autovon  587-2543,  Commercial  (315)  330-2543. 


I Introduction 


This  report  describes  the  work  done  on  developing  a pro- 
pagation model  for  fading  on  line-of-sight  (LOS)  links.  The 

propagation  model  has  been  developed  with  two  specific  purposes:  i 

first,  to  provide  DCA  with  an  LOS  model  that  could  be  incorporated 
into  its  existing  simulation  facilities;  and  second,  to  deter- 
mine the  amount  of  selective  fading  that  occurs  in  deep  fades. 

The  present  report  describes  the  general  model  development  and 
the  detailed  calculations  of  the  geometric-optic  ray  theory  of 
multipath. 

Fading  on  LOS  links  can  be  caused  from  a variety  of  physical 
mechanisms;  however  most  links  can  be  designed  to  avoid  these 
or  compensate  for  them.  The  atmosphere  duct,  on  the  other  hand, 
provides  a fading  mechanism  which  is  difficult  for  the  designer  to 
deal  with.  Because  of  their  variability  in  occurrence  and  constant 
movement,  they  become  one  of  the  major  causes  of  atmospheric  in- 
duced outages  on  well  designed  links.  We  have  thus  based  our  fad-  ^ 

ing  model  on  the  presence  of  a duct. 

In  section  II  of  this  report  we  give  a brief  review  of  the  ) 

evidence  that  ducts  can  lead  to  fading  on  LOS  links.  In  addition, 

we  show  that  when  these  fades  are  deep,  selective  fading  can  occur.  j 

In  section  III  we  develop  the  duct  model  in  detail.  The  ;| 

atmospheric  refractive  index  is  approximated  by  three  linear  segments. 

Following  this,  geometric-optics  is  employed  to  calculate  the 

i 

ray  trajectories  of  the  direct  and  multipath  rays  for  a transmitter  | 

I 
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and  receiver  having  the  same  height.  In  particular,  for  a given 
link  length,  the  number  of  contributing  multipath  rays  is  calcu- 
lated along  with  their  respective  angles  of  departure  from  the 
transmitter.  Finally  a program  has  been  developed  to  plot  the 
complete  family  of  rays  trapped  by  the  duct  to  aid  in  under- 
standing the  fading  mechanism  more  completely. 

In  section  IV  we  present  the  results  and  recommendations 


for  further  work. 


II  Channel  Modeling 

Microwave  line  of  sight  links  in  the  1 to  10  GHB  band 
suffer  from  both  slow  and  fast  fading.  The  slow  fading  has  a 
time  constant  of  several  hours  or  longer  and  is  in  the  order 
of  1 to  10  db  except  in  unusual  cases  such  as  blackout  fading. 

Past  fading,  on  the  other  hand,  is  caused  from  multipath.  Several 
different  geometric-optic  contributions  reach  the  receiver  at 
the  same  time.  One  of  them  is  the  direct  wave  contribution 
which  the  system  was  designed  to  recieve.  The  other  contributions 
to  the  received  signal  come  from  waves  that  have  been  reflected 
from  the  earth  surface  or  trapped  in  a surface  or  an  elevated  duct. 
The  time  delay  of  these  other  contributions  differs  from  the  direct 
wave  and  causes  constructive  and  destructive  interference.  If  the 
time  delay  is  substantial,  this  can  cause  deep  fast  fading  the 
received  signal.  Usually  fades  are  of  the  order  of  a second  or 
two.  In  some  cases,  fades  can  be  30  to  40  db.  deep  and  frequency 
selective. 

Although  the  fast  fading  can  occur  from  surface-reflection 
multipath  or  duct  multipath,  usually  the  surface  reflection  type 
can  be  eliminated  by  proper  design  of  the  system.  This  leaves 
duct  contributions  as  the  main  cause  of  fast  fading  in  the  micro- 
wave  frequency  range  of  interest. 

The  fact  that  multipath  was  associated  with  fast  fading  was 
documented  by  a number  of  researches  at  Bell  Telephone  Laboratories 
in  the  late  forties  and  early  fifties. In  particular,  Craw- 
ford and  Jakes ^ performed  angle  of  arrival  experiments  on  a 22 
mile  link  in  New  Jersey.  They  operated  at  about  4GH*  with 


a bandwidth  of  450  MHZ.  They  found  that  under  certain  atmospheric 
conditions,  severe  fading  occurred.  During  its  occurrence,  they 


observed  that  several  waves  were  impinging  on  the  antenna  with 

(2) 

different  angles  of  arrival.  Shortly  after  this  De  Lange  per- 
formed a short  pulse  experiment  ov?r  the  same  link  and  observed 

multipath  transmission  effects  with  path  differences  as  great  as 

(4) 

7 feet.  Much  later  in  1970  Bullington  related  the  vertical  velo- 

(5) 

city  of  the  duct  to  the  fade  rate  and  Ruthroff  demostrated  that 
application  of  simple  geometric-optic  calculations  to  a duct  could 
be  used  to  predict  Barnett's  empirical  result, i.e.  that  the  dis- 
tribution of  attenuation  is  a function  if  L^/X  where  L is  the  path 

length  and  X is  the  free  space  wavelength.  In  addition,  several 
(6,7) 

other  authors  corrobarated  the  evidence  that  in  the  absense  of 
surface  reflections,  ducts  were  the  main  cause  of  fast  fading  on 
LOS  links. 

(8.) 

In  1974,  Bello  et  al  wrote  a report  dealing  with  digital 
modem  design  for  LOS  links.  A portion  of  this  report  dealt  with 
channel  modeling.  Here  a model  for  LOS  fading  was  presented  which 
incorporated  all  of  the  physical  mechanisms  that  could  affect  LOS 
fading.  The  model  shown  in  Fig  1.,  which  is  similiar  to  the  one 
presented  by  Bello,  represents  the  transfer  function  for  a LOS 
link.  We  see  that  the  signal  leaving  the  transmitter  can  arrive 
at  the  receiver  via  the  direct,  ducted  or  surface  reflection  paths. 
Along  each  path  the  signal  is  affected  by  atmospheric  turbulence- 
labeled  as  volume  scatter  in  Fig.  1. 

The  time  varying  transfer  function  T(f,t)  relates  the  Fourier 
transform  of  the  input  and  the  output.  It  can  be  exoressed 


i 
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Multipath  Model 


analytically  as 


N -j2Tr(fo+f)Tv(t) 

T(f,t)  - £ A^(t)  [1  + g,^(t)]  e 

where  and  Tj^(t)  are  the  amplitude  and  phase  variation 

respectively  of  the  received  signals  in  the  absence  of  volume 
scattering.  They  are  usually  slowly  varying.  The  effect  of 
turbulence,  although  small,  is  added  by  the  term  gj^(t).  Its 
variation  is  usually  more  rapid.  Usually  this  term  can  be  neglected. 
The  frequency  f is  the  deviation  from  the  carrier  frequency  fo . 

Before  proceeding  with  the  detailed  calculations  of  the 
Aj^(t)  and  for  a particular  duct  we  would  like  to  demonstrate 

the  frequency  selective  character  of  the  channel  in  a deep  fade. 
First  we  assume  that  the  signal  bandwidth  is  small  enough  so  that 
exp(-j2 TrfTj^)  can  be  expanded  in  the  power  series.  Neglecting 
the  gjj(t)  for  the  time  being,  we  have 


T(f,t)  = 23 

i=0 


,(i) 


{f,t) 


(2) 


where 


,(i) 


/ \ “j2iif  (t)  2 

® (-2TTjfT3^)^ 


k=o 


n! 


(3) 


Usually  [fTj^|<<l,  so  that 

T(f,t)  = T<°)(f,t)  = 


A^(t) 


-j2TTfoTj^(t)  (4) 


We  see  in  this  case  that  the  channel  does  not  depend  oh  f and  Is 
not  frequency  selective. 

If  a deep  fade  is  encountered,  the  zero  order  multipath 
components  cancel  out  approximately  l.e., 


^f  ,t)  5:  0.  Then 


T 


^ -j2T,f,T^(t)  (5) 

k«o 


Here  we  clearly  see  that  T(f,t)  is  a linear  function  of  f and 
thus  it  is  frequency  selective.  In  addition,  we  notice  the  inter- 
esting fact,  that  since  ^ ( f ,t)  « o,  the  gjj(t)  will  now  become 
important.  We  would  thus  expect  to  see  rapid  fluctuations  in  the 
order  of  milliseconds  in  the  deep  portion  of  the  fade. 
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Ill  Multipath  Ray  Structure  for  an  LOS  Link  in  the  Presense  of  a Duct. 


In  this  section  geometric-optics  is  used  to  calculate  the 
ray  structure  of  the  wave  leaving  the  transmitter  when  a duct  is 
present.  Following  this,  the  particular  rays  that  are  intercepted 
by  the  receiver  are  identified  and  the  angle  of  departure  from  the 
transmitter  is  calculated  for  each  Intercepted  ray.  The  difficult 
task  here  is  not  the  calculation  of  the  ray  trajectories;  but  rather 
the  determination  of  which  rays  trajectories  are  intercepted  by 
the  receiver  as  a function  of  duct  height,  thickness,  receiver  and 
transmitter  height  and  length  of  link. 

A.  Problem  Formulation  - Modified  Refractive  Index 

A realistic  discussion  of  multipath  on  LOS  links  must  take 
into  account  both  the  variation  in  the  refractive  index  n(*)  where 
z is  the  height  above  the  earth's  surface  as  well  as  the  curvature 
of  the  earth.  Addressing  the  refractive  index  variation  first, 
we  assume  the  refractive  index  decreases  linearly  with  height. 

The  layer  is  incorporated  in  the  model  by  approximating  the  re- 
fractive index  by  three  linear  segments  as  is  shown  in  Fig,  2. 

Here  the  abcissa  is  measured  in  N units,  i.e. 

N(*)  - (n(2)  -1)  10* 

and  the  layer  is  located  between  h and  h+w;  the  layer  width  being  w. 
The  value  of  the  refractive  index  on  the  earth's  surface  is  taken 
as  No  and  ANi,AN2  and  AN;  are  the  absolute  values  of  the  slopes  of 
the  individual  segments.  Analytically,  N(s)  can  be  expressed  as: 


i 


P 


N,  -AN,*  , o£*<h 

N,  -AN,h  -ANtT«-h)  , h<*<h+w  (7) 

No  -AN,h-AN*w  -AN,(z-h-w)  , B ^ h+w 

Typical  values  of  these  parameters  can  be  found  in  numerous  places 

(8) 

in  the  literature  . For  purposes  of  numerical  calculations  to 
be  performed  later,  we  shall  assume  No  ■ 300  N units,  AN,  ■ 40/km 
and  160/km <AN2  < 400/km.  Note  that  ANj  must  be  greater  than  157/km 
for  ducting  or  ray  trapping  to  occur.  In  this  case  the  layer  is 
called  a duct.  Values  of  AN)  shall  not  enter  into  over  calcula- 
tions since  it  will  be  assumed  that  the  transmitter  and  receiver 
are  located  below  the  layer.  As  a result,  rays  transmitted  into 
the  third  region  never  return  to  the  reciever.  Duct  widths  are 
typically  between  40  and  300  meters. 

Before  proceeding  with  the  calculation  of  ray  trajectories, 
the  curvature  of  the  earth  must  be  taken  into  account.  This  can 
be  done  by  transforming  the  spherical  earth  with  refractive  index 
variation  n(z)  into  a flat  earth  with  a modified  refractive  index 
variation  m(z)  where 

m(z)  - n(z)  + I (8) 

with  "a"  being  the  radius  of  the  earth.  The  transformation  given 
in  (8)  is  not  exact  but  is  quite  good  for  LOS  links. 

If  eqs.  (6),  (7)  and  (8)  are  used,  we  find 


I 


i 
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m{  z) 


Wj  (z) 
(z) 


rao  + g^z  , 0£Z<h 

"*0  » h£z<h+w 


(9) 


where  the  = a”^  -AN^IO”®,  i“l,2  and  mo  = 1 + N^IO”* 

Note  that  m(z)  for  z>h+w  has  been  omitted  from  (9)  since  it  will 
not  be  used.  The  modified  refractive  index  variation  is  shown  on 
the  left  hand  side  of  Fig.  3,  It  is  plotted  in  terms  of  M units 
which  are  related  to  * by 


M = (m-l)10‘ 


(10) 


We  see  that  gi  must  be  less  than  zero  (gi<0)  for  a duct  to  exist. 
Thus  the  minimum  value  that  AN2  can  have  is  when  gj=0.  This  leads 
to  the  previously  stated  condition  that  ANj  > 157/km  for  ducting 
to  occur. 


B . Geometric-Optics 

When  the  wavelength  is  small  compared  to  the  variation  in  the 
refractive  index,  the  methods  of  geometric-optics  can  be  used  to 
obtain  approximate  solutions  to  Maxwell’s  equations.  For  LOS  links 
these  conditions  are  met.  In  addition,  since  the  waves  travel  ap- 
proximately parallel  to  the  ground,  depolarization  effects  are 
small.  It  can  be  assumed  that  each  of  the  field  components  'f  obey 
the  scalar  wave  equation: 

[7^  + k^  m^z)]')'  = 0 , k^*  w/Cq  (11) 


11- 


Ij 


where  ko  is  the  free  space  wave  number,  ui  is  the  angular  frequency 
and  Co  is  the  free  space  velocity  of  light. 

Basically,  geometric-optics  assumes  that  'i'  can  be  asymptotic- 
ally approximated  by 


iko<>(r)/  A 

'i'(r)-*  A(r)e  (l  + o(ko  ^ )I 


where  A(r)  is  the  amplitude  function  and  (r)  is  the  eikonal  or 
phase  function.  Substitution  of  (12)  into  (11)  and  comparison 
of  equal  powers  of  ko  gives  the  following  equation  for  ((i  (r)  : 


(7  'I')  ^ ■ m*  (z) 


This  equation  is  called  the  eikonal  equation.  A transport 
equation  for  A(r)  is  also  obtained. 

The  eikonal  equation  is  a nonlinear  first  order  partial 
differential  equation  that  can  be  solved  by  the  method  of  charac- 
teristics. These  characteristics  are  called  the  rays  and  they 
are  orthogonal  to  the  constant  phase  surfaces,  i.e.,  ())(£)■  const. 

For  the  special  case  of  stratified  medium,  which  we  are  con- 

(9) 

sidering,  the  ray  equation  is  particularily  simple.  It  is 


*-f  :f= 


m (z ) cosa 


(z- )-o‘ 


where  initial  coordinates  of  the  ray  and  a is 
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the  angle  between  the  x axis  and  the  ray  at  point  (x,z) , The 
function  m(z)co8a(z)  has  a constant  value  a along  any  one  ray. 
It  is  called  the  ray  parameter  and  it  serves  to  Identify  each 
ray.  The  + and  - signs  appearing  in  (14)  are  used  when  dx/dz 
(ray  slope)  is  positive  or  negative  respectively. 

Once  the  rays  have  been  found,  the  eikonal  equation  can  be 
integrated  and  the  phase  along  each  ray  can  be  calculated.  The 
result  for  the  phase  accumulation  along  a ray  from  its  initial 
point  (>t^,z^)  to  (X,z)  is 


where  ds  is  the  differential  arclength  on  the  ray.  The  time 
delay  along  this  segment  of  the  ray  is  then  simply  t = '('/Co  . 

Once  the  ray  trajectories  are  found,  they  can  be  used  in 
conjunction  with  the  transport  equation  to  find  the  amplitudes. 

The  basic  result  found  is  that  energy  Inside  a ray  tube  sur- 
rounding the  ray  is  conserved.  This  principle  can  then  be 
employed  to  find  the  amplitudes. 

The  geometric-optic  ansatz,  (12),  remains  valid  except  in 
diffraction  regions.  When  the  rays  form  a caustic  or  envelope, 
the  amplitude  A(r)  becomes  infinite  and  (12)  is  no  longer  valid. 
Such  regions  occur  for  the  problem  under  consideration,  however 
asymptotic  solutions,  valid  in  these  regions,  will  not  be  develop- 
ed here. 
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C Direct  Ray  Calculations 


We  shall  now  use  the  ray  equation  (14) presented  in  the  previous 
section  to  calculate  the  trajectory  of  the  direct  ray.  This 


is  the  ray  that  intercepts  the  receiver  whether  a duct  is  present 

or  not.  Most  LOS  links  are  designed  on  the  basis  that  this  will 

be  the  only  ray  path  from  the  transmitter  to  the  receiver.  t 

To  simplify  calculations  for  the  remainder  of  this  report 
we  will  assume  that  the  transmitter  and  receiver  are  located  at 
the  same  height.  This  height  will  be  denoted  by  Zq . In  addition,  j 

we  assume  that  they  remain  in  the  first  medium, i.e.,  Zo<h.  This 
configuration  is  shown  in  the  right  hand  side  of  Fig.  3.  The  dis- 
tance between  the  source  and  receiver  is  taken  as  L. 

The  direct  ray,  as  shown  in  Fig.  3.,  leaves  the  source  with 
a negative  slope;  it  forms  an  angle  a,  with  the  x axis.  This 

angle  can't  be  determined  now  but  will  be  found  when  L and  cxo  ■ 

are  related.  Now  employing  the  ray  equation  (14) , we  have 


An  •xamlnation  of  (17)  shows  that  as  x increases  along  the 
ray,  s decreases.  This  continues  until  the  turning  point  is 
reached.  The  turning  point,  {x^,z^) , is  the  point  at  which  a-o 
or  the  ray  has  lero  slope.  Equation  (17)  is  only  valid  for  x<x^. 
The  location  of  the  turning  point  can  be  calculated  simply  by 
using  the  constancy  of  the  ray  pareuneter  o along  the  ray.  Since 
0 must  have  the  saune  value  at  the  source  as  at  the  turning  point 

mj  (z^)cosa  ■ mj(zo)cosao  (18) 

But  a*o  at  the  turning  point,  thus 


*t  “ ("h  (zo ) ■nV))/go 


(19) 


Using  (19)  in  (17)  and  simplifying,  we  have, 


. _ i_/l  + sinaal 

‘t  To  b-BsT,"  7 


(20) 


Since  the  transmitter  and  reciever  are  at  the  same  height, 
the  direct  ray  must  be  symmetric  about  x^.  Thus  this  symmetry 
allows  us  to  determine  the  ray  trajectory  for  without 

using  (14)  again. 

In  order  to  determine  the  direct  ray  explicitly,  it  is  neces- 
sary to  ]cnow  a«  . This  can  be  obtained  from  (20)  since  L»2x^. 

We  find  . 

2jn,  (z,  )cosa  /I  + sina 

L « ml 2. 

g,  \ COSO, 


This  is  a transcendental  equation  for  , and  in  general, 
must  be  solved  numerically!  However,  in  the  case  of  LOS  links, 
a,<<l,  and  an  approximate  solution  can  be  obtained.  We  find: 


S ^ ^ (22) 

where  L/O  <<1  for  LOS  links.  The  retention  of  the  second  term 
in  the  approximation  is  necessary,  because  phase  differences  be- 
tween different  rays  are  considered,  the  effects  of  the  first 
term  in  (22)  cancel  out. 

Once  final  consideration  before  leaving  the  subject  of  the 
direct  ray  is  the  blackout  angle.  If  ao  is  too  large,  the  direct 
ray  will  be  reflected  from  the  ground  before  passing  through  its 
turning  point.  To  insure  that  this  will  not  occur,  we  require 
that  s^>o.  The  blackout  angle  can  now  be  defined  as  the  uq 
that  corresponds  to  2^*0  , i.e., 

m,{z  JcosoQg  - m,  (O)  - 

“OB  ■ 

i 

i 

We  require  that  for  the  direct  ray.  I 


(23) 

(24) 


D.  Multipath  Rays 


In  this  section  we  will  consider  the  additional  paths  along 
which  energy  can  reach  the  receiver  in  addition  to  the  direct 
ray  path.  When  a duct  is  present,  as  we  have  assumed,  some  of 
the  rays  become  trapped  in  the  duct;  for  many  receiver  locations 
several  of  these  rays  are  intercepted.  A typical  trapped  ray 
is  shown  in  Fig  3.  Our  goal  is  to  determine  the  number  of  rays 
intercepting  the  receiver  location  and  the  aj  for  each  ray.  The 
knowledge  of  ao  allows  us  to  calculate  the  phase  and  amplitude 
contributed  by  each  ray. 

Geometric  considerations  indicate  that  there  are  four  basic 
ray  types  that  contribute  to  multipath.  These  are  shown  in  Fig.  4. 
They  account  for  rays  leaving  the  transmitter  with  a positive  or 
negative  slope  and  approaching  the  receiver  with  a positive  or 
negative  slope.  Other  trapped  ray  types  consisting  of  several 
periods  of  Type  1-4  are  possible,  but  do  not  appear  to  be  as 
important  as  these. 

Before  proceeding  with  the  particulars  of  the  ray  calculations, 
we  will  find  the  conditions  for  ray  trapping  in  the  duct.  The 
duct  is  defined  as  the  region  < z < h+w  where  z^  is  defined 
from  the  equation 


mi (z^)  = m2 (h+w)  (25) 

Reference  to  the  left  hand  side  of  Fig  3 shows  that  z^  can  be 
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found  by  a simple  geometrical  construction.  To  simplify  the 
calculation  to  be  made,  we  will  assume  that  the  duct  is  always 
of  an  elevated  type,  i.e,,  Zjj>o.  This  will  eliminate  the  pos- 
sibility that  some  trapped  rays  may  be  reflected  from  the  ground. 

Next,  the  maximum  angle  of  trapping  will  be  computed.  One 
can  show  that  as  the  lag | increases,  the  upper  turning  point, 

Z.  , for  trapped  rays  increases.  Since  it  can't  be  greater  than 
h+w  we  have  j 

mi  (zo ) cos  I = m2  (h+w)  (26)  | 


where  mt  and  m2  are  defined  in  (9) . 

This  condition  came  from  the  constancy  of  the  ray  parameter  a. 
Thus 


cos  I a 


r\2  (h+w) 

mi (Zg ) 


(27) 


We  see  that  for  trapping  < ug  < Typically  is 

less  than  .5  degrees  for  standard  ducts  and  LOS  links.  An  exami- 
nation of  (25)  and  (27)  shows  that  if  = zg  then  a =0.  This 

u m 

means,  as  expected,  that  when  the  source  zg  moves  outside  the 

duct  there  are  no  rays  which  are  trapped.  In  addition,  we  can 

also  show  that  if  the  source  is  in  the  duct  ^ a • This  means  that 

OB  m 

the  direct  ray  is  not  one  of  the  trapped  rays.  This  is  an  ex- 
pected result,  since  the  direct  ray  would  remain  if  the  duct 
were  removed. 

We  will  now  focus  our  attention  on  rays  of  type  1 and  2. 

The  two  ray  types  are  mirror  imao'^s  of  each  other.  Treating 
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ray  type  1 first,  we  employ  the  ray  equation  given  in  (14)  to 
relate  ae  and  L.  We  find 


^ adz ' 

^*t2  odz* 

J /m?(z')-0*  * 

J /mf(z')-o^  ^ 

J /m?  ( z ' ) -a  ^ 

Zo 

"ti 

h 

adz-' 

J /mi(z')-a*  “ J /m? (z ’ ) -a* 

*t,  ^ 


where  mi(z^,)  ■ in2(z^^)»  mi(zo)  cos  ao  -o  with  mi  and  mj  being 
given  in  (9) . Integration  of  the  above  equations  gives 


where 


L - 


^ cosh”* (b 


secao 


). 


Oo  <o 


b 


mi  (h) 

mi  (zo ) 


(29) 


gogi 

gu+gi 


(30) 


Because  of  the  mirror  symmetry  of  type  1 and  2,  the  equation  for 
the  type  2 ray  is  the  same  as  (28)  with  the  terms  in  reverse  order 
with  ao>o.  Thus  (29)  is  the  L versus  ao  is  the  same  for  both 
ray  types. 

To  find  ao  in  terms  of  L from  (29)  we  must  again  assume  ao  is 
small.  The  approximate  expression  obtained  is 


(Typ«  1) 


a.i  - - TT  Vl*  - L*,„, 


Is  ll 


oti  - " ^IlN  <'^P* 


with 


H«IN  " g, 


ySgoAh 


Ah  ■ »#-h 


Examination  of  (31)  ahowa  that  no  raya  of  type  I or  2 contribute 
It  L < In  a later  aection,  when  the  ray  trajectories  are 

plotted  will  aae  that  the  point  L « ^ focus. 

Since  a lower  bound  to  L exists,  we  might  ask  if  there  is  a 
maximum  L.  The  answer  is  affirmative  and  it  corresponds  to  the 
ray  that  has  taken  on  the  initial  angle  |«^|>  This  maximum  dis- 
tance is  called  and  Lgj^2  types  1 and  2 respectively. 

By  using  (27)  in  (29),  we  find 

SmI  " Sm2  “ (33) 


Although  (31)  and  (33)  are  approximate  expression  for  small  angles 
of  departure  they  are  quite  accurate  for  LOS  calculations. 

The  oo  versus  L relationship  for  the  remaining  two  rays  types 
(type  3 and  4)  can  again  be  calculated  using  (14).  The  result  Isi 


L ■ cosh' ^ ^bsec a o)  - ~ cosh”  * ^seca#^  , a#>o  (341 


(Type  3) 


L ■ ^ cosh"*  (bsecoo)  ^ cosh”*  (seoou)  » ao<o  (35) 

fliji  \ ' go  ' ' 

(Type  4) 


As  before,  although  (34)  and  (35)  are  exact,  we  wish  to  find  ao 
in  terms  of  L.  Making  the  small  a*  approximation  we  find 


where 


with 


at/g. 


2t/gT 


^IN  SlIN 


' “ 03  "0 


(+) 


t - 1 


■(RJ 


We  see  that  both  type  3 and  4 rays  do  not  exist  when  Since  one 

can  show  that  t £ 1,  we  have  1 ^MIN*  ^ examination  of  (36)  shows 

that  two  rays  of  type  3 exist  if  botn  ktl  and  ko  I ax«  greater  than 
zero.  Similarllv.  two  rava  of  tvoe  4 exist  if  both  ait^  and 
are  less  than  zero.  We  find  that  one  ray  of  type  3 and  type  4 exists 
if  L > while  two  rays  of  each  type  exist  if  < L < 

llaxlmum  distances  for  rays  type  3 and  4 also  can  be  found  by 
putting  |oo|"|aj  in  (34)  and  (35).  We  find 


-BM3  ■ 
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From  the  algebraic  form  of  (33),  (40)  and  (41),  we  see  that 


Sm3  - ^M1  “ ^M2  - ^M4 

I 

The  results  of  this  section  give  us  an  explicit  method  for 
calculating  the  number  of  rays  that  contribute  at  a given  L. 

These  formula  also  give  us  the  ray  parameters  associated  with 

i 

each  contributing  ray.  Once  these  pareuneters  are  )cnown,  cal- 
I culation  of  the  time  delay  and  phase  along  each  path  follows  easily 

I from  geometric  optic  considerations  as  indicated  previously. 

t 

I E Ray  Feunilies 

; In  section  C and  D we  developed  approximate  formula  for 

i 

calculating  the  number  of  rays  contributing  to  the  received  sig- 
nal at  any  particular  value  of  L.  Although  this  is  a needed 
result,  our  picture  of  the  totality  of  rays  emitted  from  the 
source  is  incomplete.  We  need  to  have  a global  understanding 
of  how  the  rays  behave  and  %*•  need  to  determine  where  the  con- 
st Ics  and  ray  boundaries  are  located. 

i 

To  gain  this  understanding,  the  ray  equation  (14)  was  inte- 
grated and  an  IBM  370-148  computer  was  used  to  compute  the  ray 
trajectories.  These  were  then  plotted  using  a Zeta  230  plotter. 

The  results  are  shown  in  Figs  4-7.  The  plots  were  made  using 


(42)  j 

]' 
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Fig  4.  Ray  Trajectories  for  h = 105  meters 
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the  modified  refractive  index  and  thus  the  range  coordinate 
measures  distance  along  the  surface  of  the  earth. 

The  gradients  of  refractive  index  were  taken  as  AN^-40  and 
AN2*350  for  all  four  plots.  The  value  of  AN2*3S0  represents 
a fairly  steep  gradient.  However,  we  wished  to  examine  ducting 
in  a worst  case  environment  (maximum  measured  gradient  is  AN2>420) . 
The  source  was  held  at  a constant  location  of  100  meters  and  the 
duct  parzuneter  w was  kept  at  50  meters  for  all  plots.  The  only 
parameter  that  varies  from  plot  to  plot  is  the  duct  height  para- 
meter h.  The  height  increases  from  105  meter  in  Fig.  4 to  170 
meters  in  Fig  7.  Thus  the  plots  taken  together,  serve  to  show 
how  the  multipath  structure  changes  as  the  duct  rises. 

An  examination  of  the  ray  families  shows  that  each  (except 
Fig  7 ) has  two  cusps.  Each  of  these  cusps  are  formed  by  the  in- 
tersection of  two  caustics  (ray  envelopes) , The  lower  cusp  al- 
ways occurs  at  the  height  of  the  source  while  the  other  cusp 
occurs  at  a height  somewhat  above  it.  The  range  of  the  lower 
cusp  in  as  given  in  (32),  while  the  lower  caustic  associated 

with  upper  cusp  crosses  the  line  z » zq  at  given  in  (38) . 

This  is  shown  more  clearly  in  Fig  8 where  the  caustics  are  shown 
without  the  rays.  Inspecting  the  ray  feunilies  as  h Increases 
shows  that  the  upper  cusp  increases  in  height  as  h increases. 

In  addition,  both  cusps  occur  at  larger  and  larger  ranges  as  h 
increases.  This  is  the  reason  we  don't  see  the  lower  cusp  in 
Fig  7 since  it  occurs  at  such  a large  range. 

To  complete  the  interrelationship  between  section  D and  E; 
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we  have  plotted  the  extreme  trapped  raya  in  Fig  8.  These  are 
the  raya  with  the  largest  |ao  \ that  will  still  be  trapped.  As 
is  shorn  in  this  figure  the  intersection  of  these  rays  with  the 
line  coincides  with  the  outer  ray  boundary  ranges  i“l» 

2,3,4. 


IV  Conclusions  and  Recommendations 


A tractable  model  for  LOS  fading  has  been  chosen  and  the 
structure  of  multipath  rays  has  been  investigated  in  detail. 
Formula  were  determined  to  compute  the  number  of  rays  contri- 
buting at  a particular  receiver  location.  In  addition,  the 
angle  of  departure  of  each  of  these  rays  has  been  found.  This 
is  important  since  the  calculation  of  amplitudes  and  phases  or 
time  delays  becomes  a straight  forward  matter  once  the  angles 
are  known.  The  analysis  was  carried  out  for  equal  transmitter 
and  receiver  heights,  however,  the  case  of  unequal  heights 
could  be  treated  in  the  same  manner. 

To  help  understand  the  ray  structure,  families  of  rays 
were  plotted.  From  these  plots  we  see  that  multipath  rays  are 
almost  never  observed  on  links  less  than  20  kilometers  long 
while  longer  links  can  have  2 to  4 rays  contributing  in  addi- 
tion to  the  direct  ray. 

To  make  use  of  the  results  obtained  in  this  report  to  simu- 
late an  LOS  link,  the  following  tasks  would  have  to  be  performed 

1.  the  calculation  of  the  amplitudes  and  phases  would  have 
to  be  completed.  Because  the  geometric-optic  cunplitudes 
become  Infinite  at  caustics,  special  consideration  would 
have  to  be  given  to  these  regions. 

2.  Once  the  calculations  in  (1)  are  completed,  a computer 
code  for  the  transfer  function  given  in  equation  (1) 
would  have  to  be  developed. 

3.  The  dynamic  aspects  of  the  model  can  be  introduced  by 


making  the  duct  height  a slowly  varying  function  of  time. 
This  might  consist  of  a slowly  increasing  deterministic 
function  to  accoxint  for  the  gradual  rising  of  the  duct 
and  a small  slowly  varying  random  portion  to  account  for 
local  fluctuations.  Some  comparison  with  statistical 
meteorological  data  should  be  made.  It  should  be  noted 
that  for  each  new  value  of  duct  height,  the  transfer 
function  computed  in  (2)  must  be  recalculated.  The  im- 
plementation of  (2)  must  be  fast  enough  to  do  this. 

If  diversity  operation  is  to  be  simulated,  the  calcula- 
tions performed  in  this  report  will  have  to  be  extended 
to  unequal  receiver  and  transmitter  heights. 

Since  implementation  of  equation  (1)  would  require  the 
input  signal  in  the  frequency  domain,  an  FPT  would  have 
to  be  done  on  the  input  and  output  signals. 

The  effect  of  turbulence  could  be  incorporated  into  the 
model  to  determine  its  effects  when  deep  fades  occur. 
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